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ABSTRACT 

Aims. We study the physical properties and amount of molecular gas in the central regions of galaxies with active nuclei. 
Methods. Maps and measurements of the J=l-0, J=2-l, J=3-2, J=4-3 ^^CO, the J=l-0, J=2-l, and J=3-2 ^^CO 
lines in the central arcminute squared of NGC 1068, NGC 2146, NGC 3079, NGC 4826, and NGC 7469, as weU as 492 
GHz CI maps in three of these are used to model the molecular gas clouds in these galaxies. 

Results. Bright CO concentrations were detected and mapped in all five objects. In all cases, the observed lines could 
be fitted with two distinct gas components. The physical condition of the molecular gas is found to differ from galaxy 
to galaxy. Rather high kinetic temperatures of 125-150 K occur in NGC 2146 and NGC 3079. Very high densities of 
0.3-1.0 X 10^ cm"^ occur in NGC 2146, NGC 3079, and NGC 7469. The CO to H2 conversion factor X is typically an 
order of magnitude less than the 'standard' value in the Solar Neighborhood. The molecular gas is constrained within 
radii between 0.9 and 1.5 kpc from the nuclei. Within these radii, H2 masses are typically 1.2 — 2.5 x 10^ Mq. The 
exception is the (relatively nearby) merger NGC 4826 with R—0.3 kpc, and M = 3 x lO'^ M©. The H2 mass is typically 
about one per cent of the dynamical mass in the same region. 

Key words. Galaxies - individual: NGC 1068, NGC 3079, NGC 7469 - ISM - centers; Radio lines - galaxies; ISM - 
molecules, CO, C°, H2 



1. Introduction 

Molecular gas is a major constituent of the interstellar 
medium in galaxies. This is particularly true for star- 
forming complexes in the spiral arms, but strong concen- 
trations of molecular gas are also frequently found in the 
inner few kiloparsec of spiral galaxies. These concentrations 
of gas play an important role in the evolution of galaxy 
centers. They provide the material for inner galaxy star- 
bursts and the accretion of massive black holes. It is thus 
important to determine the characteristics of this molecular 
hydrogen gas (density, temperature, excitation) and espe- 
cially its amount. It cannot be observed directly, and its 
properties can only be inferred from observations of tracer 
elements, of which CO is one of the most abundant and 
easiest observable. However, the CO emitting gas is not 
in LTE, and the most commonly observed ^^CO lines are 
optically thick. We have to observe CO in various transi- 
tions to obtain reliable physical results and to break the 
temperature-density degeneracy that plagues ^^CO inten- 
sities, also in an optically thin isotope (e.g. also in ^^CO). 
The observed molecular and atomic line intensities then 
provide the essential input for further modeling. 

We have therefore observed a sample of nearby spiral 
galaxy centers in various CO transitions and in the 492 
GHz '^Pi-^Pq [CI] transition. These galaxies were selected 
to be bright at infrared wavelengths, and more specifically 
to have IRAS flux densities ii2tim > 1-0 Jy. The results 
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for seven galaxies from this sample have already been pub- 
lished. These are NGC 253 (Israel et al. 1995), NGC 7331 
(Israel k Baas 1999), NGC 6946, and M 83 = NGC 5236 
(Israel & Baas 2001 - Paper I), IC 342 and Maffei 2 (Israel 
& Baas 2003 - Paper II), M 51 = NGC 5194 (Israel et 
al. 2006 - Paper III). In this paper, we present results ob- 
tained for an additional five bright and well-studied galax- 
ies. In Table[T] we have summarized the characteristics of 
their appearance. 



At least three of the galaxies discussed in this paper 
have an active (Seyfert) nucleus (NGC 1068, NGC 3079, 
and NGC 7469), and probably also NGC 4826. Two of the 
galaxies (NGC 2146 and NGC 4826) show evidence of a 
merger event and all - including most of the galaxies pub- 
lished earher (IC 342, Maffei 2, M 83, and NGC 6946) - 
show signs of enhanced star formation activity in the recent 
past or present. NGC 4826 is as nearby as IC 342, Maffei 2, 
M 83, and NGC 6946) affording similar linear resolutions 
of 200-300 pc. The other galaxies are more distant so that 
our angular resolution of 10" — 14" corresponds to a lin- 
ear resolution no better than 800-2000 pc. Such low linear 
resolutions do not reveal much structural detail in the cen- 
tral CO distributions; millimeter array observations are re- 
quired for this. However, the multi-transition observations 
of "'^^CO and ^^CO presented here allow us to characterize 
the overall physical condition of the central molecular gas 
in ways that arer not possible otherwise. 
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Table 1. Galaxy parameters 
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Table 2. CO observations log 
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2. Sample galaxies 

2.1. NGC 1068 = Arp 37 = M 77 

This is a well-known Seyfert galaxy, relatively luminous 
at radio wavelengths, with overall and eore 1.5 GHz flux 
densities of 3.8 Jy and 1.1 Jy, respectively (Condon et al. 



1990). It has been mapped in HI by Brinks ct al. (1997, 
VLA, 8"), in the J=l-0 ^^CQ transition by Scoville et al. 
(1983; FCRAO, 45"), Young et al. (1995; FCRAO, 45"), 
Kaneko et al. (1989; Nobeyama, 17"), Heifer et al. (2003, 
BIMA, 9" X 6"), in the J=l-0 and J=2-l ^^CO transi- 
tions by Planesas et al. (1989; IRAM, 22", 13") as well 
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NGC 1068 
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Fig. 1. CO in the center of NGC 1068. Top: Observed line 
profiles. Horizontal scale is LSR velocity Vlsr in km/sec. 
Vertical scale is main-beam brightness temperature T^b in 
Kelvins, using the efficiencies listed in Table 2. Note differ- 
ent scales for ^^CO, ^^CO and [CI] respectively. Center: 
(Left) J=2-l ^^CO emission integrated over the veloc- 
ity interval 950-1350 kms^^; contours are in steps of 20 
Kkms^^; here and in all other figures, the first contour is 
equal to the step value. (Right) J=3-2 ^^CO emission inte- 
grated over the same velocity interval; contours step is 15 
Kkms"-'^. Bottom: J=3-2 '^^CO; position-velocity map in 
position angle PA = 90° integrated over a strip 20" wide 
with contours in steps of 0.1 K; east is at top. 

as in the J=2-l and J=3-2 transitions of both '^^CO and 
i^CO by Papadopoulos & Seaquist (1999; JCMT, 22", 14"). 
All maps covered more or less the same area that we have 
mapped, and show only limited detail. The presence of sig- 
nificant structure at angular sizes substantially smaller than 
these beam sizes was first suggested by Scoville et al (1983), 
using a daring deconvolution scheme. This was later con- 
firmed by the J=l-0 ^^CO millimeter array maps made at 
much higher angular resolutions by Planesas et al. (1991, 
OVRO, « 3"), Kaneko et al. (1992, NMA, « 5") and in 
particular the J=l-0 ^^CO and ^•'CO maps by Heifer & 
Blitz (1995, BIMA, ^ 4") and the J=l-0 and J=2-l ^^CO 



NGC 2146 
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Fig. 2. Center of NGC 2146, as Figure 1. Top: Observed 
line profiles. Center: (Left) J=2-l -"^^CO emission integrated 
over the velocity interval 600-1200 kms~^; contours step 
is 15 Kkms~^; (Right) J=3-2 ^^CO emission integrated 
over the same velocity interval; contour step is 20 K kms~^. 
Bottom: J=3-2 '^^CO; position- velocity map in position an- 
gle PA = 303° integrated over a strip 20" wide with con- 
tours in steps of 0.1 K; northwest is at top. 



maps by Schinnerer et al. (2000, IRAM, 0.7"). These show 
a compact (size about 5" or 400 pc) circmimuclcar source of 
molecular line emission, and bright emission extending to 
about 15" (1200 pc) from the nucleus. In lower-resolution 
maps this emission fortuitously resembles a ring, but higher 
resolutions clearly reveal CO-bright spiral arms (c/ Heifer 
& Blitz 1995). The arms dominate the ^^CO emission, but 
only in the lower J transitions. The interferometer map by 
Heifer & Blitz (1995) shows strong J=l-0 ^^CO emission 
from the arms, but nothing from the nucleus. The compact 
nucleus cannot be unambiguously identified in the J=3-2 
^■'CO map measured by Papadopoulos & Seaquist (1999) 
either. The nucleus is also dominant in the light of other 
molecules (c/ Jackson et al. 1993, Tacconi et al. 1994; Heifer 
& Blitz 1995; Usero et al. 2004) that trace high densities or 
high excitation. 
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NGC 3079 
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Fig. 3. Center o NGC 3079, as Figure 1. Top: Observed 
line profiles. Center: (Left) J=2-l ^^CO emission integrated 
over the velocity interval 950-1350 kms""'^; contour step 
is 15 Kkms^^; (Right) J— 4-3 ^^CO emission integrated 
over the same velocity interval; contour step is 20 Kkms"-^. 
Bottom: J=3-2 ^^CO; position- velocity map in position an- 
gle PA = -11° integrated over a strip 20" wide with contours 
in steps of 0.04 K; north is at top. 
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Fig. 4. Center of NGC 4826, as Figure 1. Top: Observed 
line profiles. Center: (Left) J=2-l ^^CO emission integrated 
over the velocity interval 150-650 kms~^; contour step is 
5 Kkms~^; (Right) J=3-2 ^^CO emission integrated over 
the same velocity interval; contour step is 5 kkms. Bottom: 
J=2-l ^^CO; Bottom: position- velocity map in position an- 
gle PA = 115o integrated over a strip 20" wide with con- 
tours in steps of 0.03 K; northwest is at top. 



2.2. NGC 2146 

NGC 2146 has an optically disturbed appearance, char- 
acterized by a high-surface brightness core crossed by a 
prominent dust lane, and a number of weaker extraplanar 
arms and tidal features (sec Sandage & Bedke, 1994). This 
bright core measures several kiloparsecs in diameter, and is 
a strong source of radio emission (1.4 GHz flux density «1 
Jy , see de Bruyn 1977; Condon et al. 1996; Braun et al. 
2007) against which HI is seen in absorption (Taramopoulos 
et al. 2001, Tarchi et al. 2004). Taramopoulos et al. (2001) 
have presented a VLA HI map including short spacings 
that shows a neutral gas filament extending from the main 
body of galactic HI to the south over half a degree (160 
kpc) , as well as a large amount of gas extending out of the 
galaxy to the north. The HI mass of the central region (in- 
ner few arcmin) is about 1.6 x 10^ M©. Taramopoulos et 



al. (2001) conclude that NGC 2146 suffered a merger or 
close encounter in which the intruder was shredded, about 
800 million years ago. The bright core is the site; of a pow- 
erful starburst. Radio observations reveal numerous super- 
nova remnants, HII regions and water masers throughout 
(Kronberg & Biermann, 1981; Zhao et al. 1996; Lisenfeld et 
al. 1996; Tarchi et al. 2000, 2002); a low-luminosity AGN 
may be present (Inui et al. 2005). X-ray observations also 
reveal a starburst-driven superwind carrying material away 
from the galaxy's plane (Armus et al. 1995; Delia Ceca et 
al. 1999). It has been mapped in the J=l-0 ^^CO line by 
Young et al. (1988b, FCRAO, 45"; OVRO, 7") and .lackson 
& Ho (1988, Hat Creek, 6 x 5"), in the J=l-0 and J=2- 
1 transitions of both ^^CO and ^^CO by Paglione et al 
(2001, FCRAO, 45") and Xie et al (1994, FCRAO, 23") 
respectively, and in the ,7=3-2 ^^CO line by Dmnke et al. 
(2001, HHT, 24"). The highest resolution observations have 
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Table 3. ^^CO and [CI] observations log 
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Fig. 5. Center of NGC 7469, as Figure 1. Top: Observed 
line profiles. Center: (Left) J=2-l ^^CO emission integrated 

over the velocity interval 4650-5050 kms"-'^; contour step is 
5 Kkms^^; (Right) J=3-2 ^^CO emission integrated over 
the same velocity interval; contour step is 7.5 Kkms"-^. 
Bottom: J=3-2 ^^CO; position- velocity map in position an- 
gle PA = 128° integrated over a strip 20" wide with con- 
tours in steps of 0.025 K; east is at top. 



been published by Greve et al. (2006) using both the IRAM 
30 m (J=l-0, J=2-l i^cO; 22" and 13") and PdB array 
telescopes (J=l-0 i^CO and ^^CO; J=2-l ^^CO; 4-5" and 



3"). 



2.3. NGC 3079 

This is a large and bright spiral seen almost edge-on. Its 
central region is a strong source of radio continuum emis- 
sion (de Bruyn, 1977). It has been mapped in the radio at 
various resolutions - sec references in Condon et al. (1990). 
These authors list 1.4 GHz flux densities of 0.84 Jy and 
0.55 Jy for the whole galaxy and the central region respec- 
tively (see also Irwin & Saikia 2003). The center contains a 
nuclear source (AGN surrounded by water masers, Trotter 
et al. 1998; Hagiwara et al. 2004) and two wind-driven out- 
flow bubbles (Veilleux et al. 1994; Cecil et al. 2001); it is 
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a significant X-ray source (Pietsch et al. 1998). The cen- 
tral region also hosts a powerful starburst. In the plane 
of the galaxy, the cavity surrounding the AGN is marked 
by near-infrared line emission from shocked H2 extending 
out to 300 pc (Israel et al. 1998; see also Hawarden et al. 
1994). At the cavity wall, hot dust reaches evaporation tem- 
peratures (Israel et al. 1998, see also Armus et al. 1994). 
Molecular gas has been mapped with a single dish in the 
J=l-0 line by Tinncy et al. (1990. NRAO, 55") and 

Young et al. (1995, FRAO, 45"), in the J=l-0 and J=2-l 
i^CO lines by Braine et al. (1997), in the J=l-0 ^^CO and 
i^CO hues by Paglione et al. (2001, FCRAO 45") and with 
mm arrays in the J=l-0 ^^CO transition by Young et al. 
(1988a, OVRO, 8" x 6"), Sofue & Irwin (1992, NMA, 4"), 
Tacconi et al. (1996, PdB, 2"), Sofue et al. (2000, NMA, 
1.5"); Koda et al. (2002, NMA, 1.5"). HI maps have been 
presented by Braun et al. (2007). [CI] emission was mea- 
sured by Gcrin & Phillips (2000), and Israel & Baas (2002). 
Absorption has been mapped in various lines (HI: Pedlar 
et al. 1996; OH: Hagiwara et al. 2004; CH3OH: Impelhzzeri 
et al. 2008). All data indicate the presence of dense molec- 
ular material extending out to almost a kiloparsec from the 
nucleus. 
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2.4. NGC4826 = M 64 



This is a small and nearby early-type (Sab) galaxy. It is 
relatively isolated although it seems to be part of the loose 
CVn I group. Its remarkable appearance, with a very promi- 
nent and wide dust lane covering the northern edge of its 
bulge, has caused it to become known as the 'evil-eye' 
galaxy. The discovery (Braun et al. 1992) that its inner 
parts rotate counter to its outer parts, i.e. that its rota- 
tion reverses sign at about 30" (600 pc) from the nucleus, 
has led to considerable interest in this galaxy since. This 
particular kinematical signature, subsequently studied by 
Rubin (1994), Walterbos et al. (1994) and Rix et al. (1995), 
is now generally assumed to mark the acquisition of ex- 
ternal material in a merger. It is the outer gas disk that 
'counter-rotates' (Rix et al. 1995), and the inner gas disk 
may have contracted from gas that was originally at greater 
radii but lost much of its angular momentum in the merger 
(Sil'chenko 1996). The central region is a weak source of 
radio emission {Si,5ghz — 100 mJy). HI maps have been 
presented by Braun et al. (2007) and Haan et al. (2008) who 
found a total HI mass Mm = 1.8 x 10*^ Mq. The galaxy 
was mapped in ^^CO and observed in ^'^CO in both the 
J=l-0 and th J=2-l transition by Casoli & Gerin (1993), 
further studied in HI and in J=3-2 ^^CO by Braun et al. 
(1994). High-resolution millimeter array observations have 
been provided in the J=l-0 ^^CO line by Sakamoto et al. 
(1999, NMA, « 4"), in the J=l-0 and J=2-l ^^CO lines 
by Garci'a-Burillo et al (2003, PdB, 1-3") and in the J==l- 
i^co i^CO and HCN lines by Heifer et al. (2003) and 
Rosolowsky & Blitz (2005). 



2.5. NGC 7469 = Arp 298 



This is the most distant galaxy in our sample. It is a classi- 
cal Seyfert 1 galaxy, and it is well-studied at various wave- 
lengths. It appears to be interacting with its companion 
IC 5283, 1.3' to the NNE. It should be noted that vari- 
ous distances quoted in the literature are about twice the 
distance assumed in this paper. High-resolution VLA radio 
maps have been presented by e,g. Ulvestad et al. (1981); 
Condon et al. (1982, 1990) and Wilson et al. (1991). The 
radio emission (total flux density Si,^ghz — 160 mJy) ex- 
tends over the central k, 8" (1150 pc) and arises mostly 
from a circumnuclear starburst 'ring' (Wilson et al. 1986, 
1991) surroimding a compact nuclear radio source (flux 
density S^^ghz = 21 mJy). The nucleus is a variable X-ray 
source (Barr, 1986; Fabbiano et al. 1992). The circumnu- 
clear starburst dominates the energy budget as it represents 
two thirds of the bolometric luminosity of the entire galaxy 
(Genzel et al 1995) and is on the order of !0^^ Lq. The 
galaxy is a powerful object at infrared wavelengths and was 
mapped with millimeter arrays in the J=l-0 ^^CO line by 
Sanders et al. (1988, OVRO, 6") and Meixner et al. (1990, 
Hat Creek, 2"), as well as more recently in the J— 2-1 ^^CO 
and J=l-0 HCN hues by Davies et al. (2004, PdB, 0.7" and 
2" respectively). Significant detail is, however, only revealed 
by the very high-resolution J=2-l ^^CO map. 



3. Observations and results 

3.1. Observations 

The observations described in this paper were carried out 
with the 15m James Clerk Maxwell Telescope (JCMT) on 
Mauna Kea (Hawaii) Q and with the IRAM 30m telescope 
on Pico Veleta (Spain). 

3.1.1. JCMT Observations 

At the epoch of the mapping observations (1991-2001) the 
absolute pointing of the telescope was good to about 3" 
r.m.s. as provided by pointing observations with the JCMT 
submillimeter bolometer. The spectra were calibrated in 
units of antenna temperature T^, correcting for sideband 
gains, atmospheric emission in both sidebands and tele- 
scope efficiency. Calibration was regularly checked by ob- 
servation of a standard line source. Further observational 
details are given in Tables[2]and[31 Most of the observations 
were carried out with the now defunct receivers A2, B3i 
and C2. Observations in 2001 were obtained with the cur- 
rent receivers B3 (330/345 GHz) and W/C (461 GHz). FuU 
details on these receivers can be found at the JCMT web- 
site (http://docs.jach.hawau.edu/JCMT/HET/GUIDE/). 
Up to 1993, we used a 2048-channel AOS backend cover- 
ing a band of 500 MHz (650 km s'^ at 230 GHz). After 
that year, the DAS digital autocorrelator system was used 
in bands of 500 and 750 MHz. Integration times (on+off) 
given in Tables[2]are typical values appropriate to the maps. 
We subtracted second order baselines from the profiles. All 
spectra were scaled to a main-beam brightness tempera- 
ture, Tnib = Tp^/rj^To] values for the values of Ty^b used are 
given in Tables [^151 

All maps were made in rectangular grids (parameters 
given in Table[2), rotated to parallel the major axis of the 
galaxy mapped. The maps shown in this paper all have been 
rotated back to a regular grid in right ascension and decli- 
nation. Some additional remarks are in order. In the case of 
NGC 1068, we combined our data as listed in Tabel[2]with 
those obtained by Papadopoulos & Seaquist (1999), taken 
from the JCMT archives, and re-reduced. The CO maps 
presented in this paper are based on both datasets. The ve- 
locity widths of the lines from the edge-on galaxy NGC 3079 
were similar to the bandwidth of the JCMT backends used. 
For this reason, we observed the lines at three slightly differ- 
ent velocity settings. The resulting spectra were combined 
and concatenated to yield the final spectra. 

3.1.2. IRAM observations 

Observations of both ^^CO and ^^CO in the J=l-0 and 
J=2-l transitions at resolutions of 21" and 12" respec- 
tively were obtained towards the centers of all galaxies with 
the IRAM 30m telescope in the period 2005-2007. We used 
the IRAM facility low-noise receivers A and B in both the 
3mm and the 1mm bands. For backends we used the IMHz 



The James Clerk Maxwell Telescope is operated on a 
joint basis between the United Kingdom Particle Physics and 
Astrophysics Council (PPARC), the Netherlands Organisation 
for Scientific Research (NWO) and the National Research 
Council of Canada (NRC). 

^ The IRAM 30m telescope is supported by INSU/CNRS 
(France), MPG (Germany), and IGN (Spain). 
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Table 4. Measured central ^^CO and ^^CO line intensities 
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Notes:" beams not directly observed but obtained from convolved maps (see Section 3.2) are in italics; measured position is 10" 
off nucleus. 



and 4MHz filter banks, as well as the VESPA correlator. 
Initially, all four transitions were observed simultaneously, 
therefore with identical pointings. After this, the signal-to- 
noise ratio of the ^^CO profiles was increased by additional 
simultaneous observation of the two transitions, where we 
took care to ensure that the pointing was the same as for 
the earlier observations. Calibration and reduction proce- 
dures were standard and similar to those described in the 
previous section. The assumed jymb values are likewise given 
in Tables[2] and [S] The IRAM observations were especially 
important not only because they provide measurements of 
the 3mm J— 1-0 transitions not possible with the JCMT. 
The IRAM aperture is twice that of the JCMT, so that the 
combination of the two datasets provides J—2-1/ J=l-0 and 
J=3-2/J—2-l ratios in closely-matched beams, 

3.2. Results 

Central line profiles in all observed transitions are shown 
in Figs.[T] through O We mapped the distribution of the 
molecular gas in the galaxy centers in the J=2-l, J=3- 
2, and J=4-3 transitions, typically over an area of about 
one arcminute across. These maps are also shown in the 
Figures, as are the major axis position velocity diagrams of 
the central CO emission. 

Measured line intensities towards the center of each 
galaxy are listed in Tabled) This Table contains entries both 
at full resolution, and at resolutions corresponding to those 
of transitions observed in larger beams. These entries were 
extracted from maps convolved to the resolution quoted. 
The convolution was carried out with the map interpola- 
tion function of the SPECX data reduction package^. This 
function interpolates, where necessary, the map intensities 
by using adjacent pixels out to a preset distance (usually 2 
beams) and convolves the observed and interpolated data 
points with a two-dimensional Gaussian function ('beam') 
of predetermined halfwidth. This halfwidth is chosen such 



http://docs.jach.hawaii.edu/JCMT/cs/005/ll/Iitml/nodel29 



that the actual observing beam convolved with the 2D- 
Gaussian would yield a beam corresponding to the desired 
resolution. Velocities in a position-velocity map can also be 
'smoothed' by applying a convolving Gaussian function in 
a similar manner. 

Finally, we have determined the line intensity ratios of 
the observed transitions at the center of the galaxies ob- 
served. These line ratios have been corrected for small off- 
sets between the various datasets. They have been deter- 
mined independently from the intensities given in Tabled) 
Specifically: (a) we used data taken at the same observing 
run wherever possible, (b) we determined ratios at vari- 
ous resolutions convolving the maps to the relevant beam 
widths, (c) we verified ratios by comparing profile shapes 
in addition to integrated intensities, and (d) the line ratios 
have been corrected for small offsets between the various 
datasets. In Table 5, whenever data were insufficient for 
convolution, or altogether lacking (e.g. [CI]), no entry is 
provided. 



4. Analysis 

4.1. Radiative transfer modeling of CO 

We have modeled the observed ^^CO and ^'^CO line 
intensities and ratios with the large velocity gra- 
dient (LVG) radiative transfer models described by 
Jansen (1995) and Jansen et al. (1994) - but see 
also Ho gerheijde fc van der Tak ( 2000) and the 
web page ]http://www. strw.leidenuniv.nl/ 'michiel/ratran/. 
These codes provide model line intensities as a function of 
three input parameters per molecular gas component: gas 
kinetic temperature Tk, molecular hydrogen density n{H2) 
and the CO column density per unit velocity iV(CO)/dy. 
By comparing model to observed line ratios, we may iden- 
tify the physical parameters best describing the actual con- 
ditions at the observed positions. Beam-averaged properties 
are determined by comparing observed and model intensi- 
]hiaal, In principle, with seven measured line intensities of two 
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Table 5. Adopted velocity- integrated line ratios in galaxy centers 
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Table 6. Radiative transfer model parameters 
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isotopes, properties of a single gas component are overde- 
termined as only five independent observables are required. 
In practice, this is mitigated by degeneracies such as oc- 
cur for ^^CO. In any case, we find that fits based on a 
single-component gas arc almost always incapable of fully 
matching the data. However, we usually obtain good fits 
based on two gas components. The solutions for a two- 
component gas are slightly underdetermined but we can 
succesfully compensate for this by introducing additional 
constraints as described below. The physical gas almost 
certainly has a much wider range of temperatures and den- 
sities. However, with the present data two components is 
the maximum that can be considered, as adding even one 
more component increases the number of free parameters to 
the point where solutions arc physically meaningless. Thus, 
our two-component model gas is only an approximation of 
reality but it is decidedly superior to one-component model 



solutions. As long as significant fractions of the total molec- 
ular gas mass are limited to similar segments of parameter 
space, they also provide a reasonably realistic model for the 
actual state of affairs. Specifically, our analysis is less sen- 
sitive to the occurrence of gas at very high densities and 
temperatures, but such gas is unlikely to contribute singif- 
icantly to the total mass. 

In order to reduce the number of free parameters, we 
assumed identical CO isotopical abundances for both gas 
components. Furthermore, in a small number of starburst 
galaxy centers (NGC 253, NGC 4945, M 82, IC 342, He 2- 
10), values of 40 ± 10 have been suggested for the isotopical 
abundance [i2C0]/[i^^C0] (Mauersberger & Henkel 1993; 
Henkel et al. 1993, 1994, 1998; Bayet et al. 2004), some- 
what higher than the characteristic value of 20-25 for the 
Milky Way nuclear region (Wilson & Rood 1994). We there- 
fore have adopted an abundance value of [^^CO]/[^'^CO] = 
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Table 7. Beam-averaged model parameters 
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40 in our models. We identified acceptable fits by search- 
ing a grid of model parameter combinations (10 K < Tk <^ 
150 K, lO^cm^^ < „(H2) < lO^cm-^, 6 x lO^^cm-s < 

N{CO)/dV < 3 X 10^^ cm"^) for model line ratios match- 
ing the observed set, with the relative contribution of the 
two components as a free parameter. Solutions obtained in 
this way are not unique, but rather define a limited range of 
values in distinct regions of parameter space. For instance, 
variations in input parameters may to some extent compen- 
sate one another, producing identical line ratios for some- 
what different combinations of input parameters. Among 
all possible solution sets, we have rejected those in which 
the denser gas component is also much hotter than the more 
tenuous component, because we consider the large pressure 
imbalances implied by such solutions physically implausi- 
ble, certainly on the kiloparsec scales observed. 

The results of our model fitting procedure are summa- 
rized in Table[6] 

4.2. Beam-averaged molecular gas properties 

The chemical models presented by van Dishoeck & Black 
(1988) show a strong dependence of the iV(C)/A^(CO) col- 
umn density ratio on the total carbon (C-I-CO) and molecu- 
lar hydrogen column densities. If line intensities of all three 
species, CO, [CI] and [CII] are known, minimal assump- 
tions suffice to deduce total carbon and hydrogen column 
densities and masses. 

If C° and C+ have not been observed, we can still re- 
late the radiative transfer model results to molecular hydro- 
gen column densities by taking the [C]/[H] gas-phase abun- 
dance ratio as an additional constraint. We must then find 
the molecular hydrogen column density that satisfies both 
the total carbon iVc = N{C) + N{CO) column-density 
and the carbon ratio fc — N{C)/N{CO) commensurate 
with the above-mentioned models and [C]/[H] abundance. 
We intend to compare the results for various galaxies, sev- 
eral of which are lacking either [CI] or [CII] measurements, 
or both. Thus, for the sake of consistency, we use this sec- 
ond method for all galaxies, whether or not [CI] or [CII] 
have been detected. 

Oxygen abundances and gradients of many nearby 
galaxies have been compiled by Vila-Costas & Edmunds 
(1992), Zaritsky et al. (1994), and Pilyugin, Vflchez & 
Contini (2004). Extrapolated to zero radius, the first two 
compilations imply nuclear metalicities many times solar. 
However, Pilyugin et al.'s more recent compilation yields 
nuclear metalicities on average only 1.5 times solar. These 



are nevertheless quite uncertain, because the galaxy disk 
values and gradients were derived from HII regions at vari- 
ous distances to the center, and it is not very clear whether 
such results can be extrapolated to the rather different 
galactic center conditions. Another practical complication 
is that of our sample, only NGC 1068 is included in these 
compilations, although we have an estimate for NGC 7469 
from X-ray observations (Blustin et al. 2007). In the fol- 
lowing, we have conservatively assumed a metalicity twice 
solar to apply to the galaxy central regions. The results 
published by Garnett et al. (1999) suggest very similar car- 
bon and oxygen abundances at these metalicities, so that 
we adopt [C]/[H] w 1.0 x 10^'^. As a significant fraction 
of carbon is tied up in dust particles and thus unavailable 
in the gas-phase, we have also adopted a fractional cor- 
rection factor 6c = 0.27 (see for instance van Dishoeck & 
Black 1988), so that iVn/A^c = [2N{ B.2) + N{m)]/[N{CO) 
+ iV(CII) + iV(CI)] 3700, uncertain by about a factor of 
two. In Table[7]we present beam-averaged column densities 
for CO and C (=C°-|-C+), as well as H2 derived under the 
assumptions just discussed. As the observed peak CO in- 
tensities are significantly below the model peak intensities, 
only a small fraction of the (large) beam surface area can 
be filled with emitting material. 

Although the analysis in terms of two gas components is 
superior to that assuming a single component, it is still not 
fully realistic. For instance, the assumption of e.g. identical 
beam filling factors for the various species ( ^^CO, ^'^CO, 
C°, and C"*") is not a priori plausible. Fortunately, these 
assumptions are useful but not critical in the determina- 
tion of beam-averaged parameters. If, by way of example, 
we assume a smaller beam filling factor, the model cloud 
intensity increases. This generally implies a higher model 
column-density which, however, is more strongly diluted 
by the beam. The beam-averaged column-density is mod- 
ified only by the degree of non-linearity in the response 
of the model parameters to a change in filling factor, not 
by the magnitude of that change. A similar state of af- 
fairs exists with respect to model degeneracy, i.e different 
model parameter combinations yielding (almost) identical 
line ratios. Although, for instance, the model fractions of 
hot or dense gas may differ, the final derived beam-averaged 
molecular hydrogen densities are practically the same. In 
fact, the assumed gas-phase carbon abundance dominates 
the results. The derived beam-averaged column-densities 
(and projected mass densities) roughly scale inversely pro- 
portional to the square root of the abundance assumed: 
Nav oc [C]/[H]gas~'^'^ ■ As our estimate of the central ele- 
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mental carbon abundance, and our estimate of the carbon 
depletion factor are unlikely to be wrong by more than a 
factor of two, the resulting uncertainty in the final values 
listed in Table[7] should be no more than 50%. 

5. Discussion 

5.1. NGC 1068 

The extended inner star-forming zone has line ratios that do 
not allow for much variation of the gas properties. Given 
the size (deconvolved dimensions of 27 x 24") of the re- 
gion, its general appearance, and the relatively slight con- 
tribution of the compact circumnuclear disk, the emission 
in the lower CO transitions refers to the extended emis- 
sion associated with the star-forming spiral arms and disk 
(see, for instance, Wynn- Williams, Becklin & Scoville 1985 
or Le Floc'h et al. 2001). The major axis- velocity map 
shown in in Fig.[T] illustrates the dominance of the spiral 
arms in the J=3-2 ^^CO distribution in the shape of the 
bright maxima at the extreme velocities, with much less 
bright CO emission closer to the nucleus. The central min- 
imum is not pronounced, however, because it is filled in 
by emission from the unresolved circumnuclear source (see 
also Fig. 6 in Heifer & Blitz 1995). We did not use the 
J=4-3 results in our analysis. The models predict a J=4- 
3/J=2-l line ratio of about 0.3, only half of that actu- 
ally observed fTable [3?2)) . However, much of the emission 
in this transition arises in the compact circumnuclear gas, 
rather than in the extended source (see Fig.[6l We find a 
relatively cool molecular gas (Tkin ~ 30 K) with densi- 
ties ranging over (riHa = lO'^-lO'* cm^'^). The inner star- 
forming zone shows up quite clearly in the high-resolution 
HI map by Brinks et al. (1997), reaching peak column den- 
sities N{HI) « 3 X 10^1 cm-2 in 8" beams. Such H2 and 
HI values are typical values for photon-dominated regions 
(PDRs) associated with star formation. 

Most, but not all (60%) carbon must be locked up in 
CO. The overall filling factor of the molecular gas is 7%. 
The molecular gas mass mapped over the central 2R = 2 
kpc is slightly more than M-^^ = 10* Mq, i.e a bit less 
than one per cent of the dynamical mass M « 1.5 x 10^° 
implied by the major-axis- velocity map in Fig[TJ 

The compact circumnuclear gas is reasonably well re- 
solved by the 10"-14" beam data, but the lack of ^^CO 
intensities severely hamper attempts to model it. The ra- 
tio i^cO/^^CO > 6 implied by the BIMA observations 
by Heifer & Blitz (1995) does not present a significant 
constraint. Nevertheless, the models do imply that a sig- 
nificant fraction (10-50%) of this circumnuclear gas must 
have a density n,^^ w 10^ cm~^ and a temperature of 
30-60 K. This is consistent with the values estimated by 
Jackson et al. (1993), Tacconi et al. (1994), Heifer & Blitz 
(1995), Papadopoulos & Seaquist (1999), and Usero et al. 
(2004). The remainder of the gas must be rather tenuous 
(riHa = 100 — 1000 cm^'^) at an undetermined tempera- 
ture. The emission from this gas is particularly prominent 
in the J=4-3 ^^CO and [CI] maps shown in Fig. [SI and the 
maps by Jackson et al. (1993) and Heifer & Blitz (1995) 
show a similar prominent role for J=l-0 HCN. The lack 
of "'^^CO emission precisely at the nucleus where the [CI] 
emission is strongest, implies that the actual [CI]/ ^^^00(2- 
1) ratio is significantly higher than the beam-confused value 
given by Israel & Baas (2002) - already one of the highest 




Fig. 6. Central 30" of NGC 1068. Top left J = 3 - 2 ^^CO 
with contour step of 10 K kms~^, top right J = 3 — 2 ^■^CO 
with contour step of 1 Kkms^^, bottom left J = 4 — 3 
^■^CO with contour step of 20 Kkms~^, bottom right [CI] 
with contour step of 10 Kkms~^. Note dominance of the 
unresolved nuclear source in the latter two maps 

in their galaxy sample. Inside the star-forming spiral arms, 
and thus also towards the nucleus, there is very little HI 
{N{HI) < 2 X 10^0 cm-2 cf Brinks et al. 1997). Much of 
the circumnuclear gas thus is dense and warm, and appears 
to be excited by the AGN rather than by luminous stars 
(c/Lepp & Dalgarno 1996; Usero et al. 2004, Meijerink et 
al. 2006, 2007). Strong neutral carbon emission is charac- 
teristic for X-ray excitation (Meijerink et al. 2007). 

The extended bar. Particularly interesting is the struc- 
ture exhibited by [CI]. The unresolved intense source, coin- 
ciding with the nucleus, is accompanied by tongues of [CI] 
emission stretching out to the northeast and southwest in 
a position angle of about 45°. The two extensions are real, 
and not due to instrumental effects, as the [CI] profiles in 
the map accurately represent the rapidly changing veloc- 
ity field. This extended [CI] emission is in the same posi- 
tion angle as the northeastern radio jet (Wilson & Ulvestad 
1987), and coincides with the northeastern near- infrared 
[Fell] line and X-ray continuum emission extensions (Blietz 
et al.l994; Young et al. 2001). However, the carbon emission 
in our map also extends in the opposite direction and thus 
more closely corresponds to the distribution of the near- 
infrared K-band continuum emission (Rotaciuc et al. 1991; 
Davies, Sugai & Ward 1998), the near-infrared polarized 
emission (Lumsden et al. 1999), and the J=l-0 ^^CO emis- 
sion in the fuUy-UV-sampled array map by Heifer & Blitz 
(1995). Interestingly, the southwestern molecular gas exten- 
sion clearly seen in their map is absent in the ^^CO maps 
by Schinnerer et al. (2000) and Usero et al. (2004) that 
lack the short-spacing information. There is no HI emission 
corresponding to the [CI] emission (c/ Brinks et al. 1997), 
although some HI absorption {N{HI) —l — 2x 10^^ cm~^ 
in a w 2" beam) is seen against the continuum of the south- 
western bar. The nature of the extended carbon emission is 
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Fig. 7. Central 30" of NGC 2146. On the left J=4-3 ^^CO, 
and on the right J— 6-5 ^^CO. Both maps have contour 
steps of 10 K kms"^. The relatively high intensity of J=6-5 
"'^^CO confirms that much of the molecular gas in NGC 2146 
is at elevated temperatures. 

not clear. Its presence at directions devoid of X-ray or [Fell] 
emission suggests that the observed strong [CI] emission 
does not originate in an X-ray-excited dense gas. Rather it 
appears that this emission arises from (photo) dissociated 
gas streaming along the stellar and molecular bar discussed 
by Tacconi et al. (1994) and Helfer& Blitz (1995). Its na- 
ture would thus be different from that of the neutral carbon 
in the circumnuclear source. 

5.2. NGC 2146 

Our maps show bright CO extending over 37" ( J=3-2 size) 
coincident with the bright central starburst region. The de- 
convolved extent of 28" (2.4 kpc) is identical to that found 
by Dumke et al. (2001) with a larger beam, and indeed to 
the extent of the resolved ^^CO and dust emission observed 
at much higher resolution with the various millimeter ar- 
rays (Young et al. 1988b; Jackson & Ho 1988: Greve et 
al. 2006). Thus, there appears to be very little molecular 
gas beyond radii of a kiloparsec. The CO extent also corre- 
sponds quite closely to that of the bright radio emission in 
the VLA maps by Kronberg & Biermann (1981) and Zhao 
et al (1996); the array of discrete radio sources mapped at 
very high resolution (0.4") by Tarchi et al. (2000) covers 
a full span of 25". The appearance of the high-resolution 
CO maps, the unmistakable double structure in maps of 
peak rather than integrated brightness temperature (not 
shown) and the shape of the CO distribution in the major 
axis- velocity diagram (Fig.[2|), exhibiting minimum inten- 
sity at the galaxy center, establish that the observed CO 
does not homogeneously fill the central volume. Rather, it 
must be distributed in a warped ring (Greve et al. 2006), or 
much more likely, in CO-enhanced spiral arms seen mostly 
edge-on. The central few hundred parsec should be largely 
free of gas. The ionized gas measured by Zhao et al. (1996) 
appears to cover the full range of velocities seen in Fig. [2] 
within the central 5". Tarchi et al. (2002) found two water 
kilomasers marking star- forming regions. One is located in 
the northwestern CO emitting region, 7 from the center 
with Vlsr ~ 825 kms~^, and the other projected very 
close to the nucleus with Vlsr = 1010 kms~^, also within 
the CO contours in Fig[2]-bottom. 

Our LVG analysis effectively yields two different model 
solutions (see Table [6]). The first one combines quite warm 
{Tkin = 125 K) and not very dense {n-a^ = 10^ cm~^) 
molecular gas with a three times smaller amount of cooler 
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Fig. 8. Major-axis-velocity map of the [CI] emission from 
NGC 3079. Contour step is 40 mK. Note relative minimum 
of [CI] emission in the center of NGC 3079. 



{Tkin — 30 K), much denser (riHs — 10^ cm^"^) gas. The 
second solution finds essentially all molecular gas to be 
warm {T^m — 150 K), but most of it is not very dense 
(riHa — 10"^ cm~'^) with a small fraction reaching den- 
sities more than an order of magnitude higher (nH2 = 
3 X 10" cm-3). 

However, we have also detected J=6-5 ^^CO emission 
from the center of NGC 2146, at intensities of about 0.5 
to 0.7 times that of the J=4-3 ^^CO emission in the same 
beam (Figure[7]). The second of the two model solutions 
('all gas is warm') is in much better agreement with this 
ratio than the first. In terms of overall results. As Table[7] 
shows, the difference is not great in terms of beam-averaged 
results, as indeed should be expected from the discussion 
in Section 4. 

5.3. NGC 3079 

Our ^^CO maps (see Fig.|31) show a just-resolved CO core 
superposed on the more extended emission of the edge- 
on galaxy. In Fig[3]-bottom the rapidly rotating gas clearly 
stands out from the more sedately moving gas in the disk 
of NGC 3079. The central velocity field is complex. The 
dynamics of the central molecular gas of NGC 3079 have 
been discussed in considerable detail by Sofue et al. (2001) 
and Koda et al. (2002) who propose the presence of both a 
weak bar, and an extremely massive black hole in the cen- 
ter. The diameter of the molecular core in the galaxy plane 
is consistent with the diameter of 14" derived from millime- 
ter array maps (e.g. Sofue et al. 2001). In the J=4-3 ^^CO 
line, the core appears to be more compact than in the lower 
transitions, but this just reflects the contribution from the 
large-scale disk weakening with increasing transition. 

The behavior of the neutral carbon emission differs from 
that in NGC 1068 described above. In NGC 3079, the [CI] 
emission does not peak at the center of the galaxy. Rather, 
as Fig. [8] shows, it peaks at either side of the galaxy, and 
the emission from the center is relatively weak. This sug- 
gests a distribution in which a larger fraction of the central 
region is free of neutral carbon than of carbon monoxide. 
As was the case for NGC 2146, the observed line ratios al- 
low two essentially different LVG model solutions. The first 
one provides a rather good fit to the data and puts a very 
large fraction of the molecular gas in a hot {Tkin = 150 
K), tenuous (71H2 = 10^ cm~^) gas, leaving a relatively 
small fraction of the gas in a cool {Tkin = 30 K) and very 
dense (nna = 10^ cm~^) molecular gas. The second so- 
lution provides a fit that is not as good, but still within 
the line ratio uncertainties. Now, slightly less than half of 
the molecular gas is hot {Tkin — 150 K) and modestly 
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dense {nn^ — 10'^ cm^'^), whereas the remainder is cold 
(Tkin = 20 K) and reasonably dense (nH2 = 10'' cm~'^). 
Braine et al. (1997) found a mean dust temperature of 30 
K, which suggests that the second model solution is more 
likely than the first. The mean molecular gas column den- 
sities of NGC 3079 are highest of all the sample galax- 
ies, which is not surprising in view of the almost exactly 
edge-on orientation of the galaxy, and the high extinction 
(Ay > 6™, Israel et al. 1998) towards the very center. The 
central molecular mass is « 2 x 10* Mq, for either 

model solution. This mass is much lower than the values de- 
rived by the various authors quoted above using 'standard' 
CO-to- H2 conversion factors. It is, however, very close to 
the mass derived by Braine et al. (1997) using a different 
method. An important consequence of this lower molecular 
gas mass determination is that the gas no longer constitutes 
such a large fraction (half!) of the dynamical mass. 

5.4. NGC 4826 

The maps in Fig.[3]show that essentially all molecular gas is 
contained with a distance of 30" (600 pc) from the nucleus 
- i.e. all molecular gas is located within the radius at which 
the rotation curve reverses sign. This is particularly well 
illustrated by the J=2-l ^^CO major-axis velocity map at 
bottom. Unlike the p-V maps of the other galaxies, this 
map shows very clear boundaries at either side. Although 
the velocity-integrated CO emission appears to be peaked 
on the center in all maps, the position-velocity maps in 
Figs.[3]and[n]show that this is not the case. In all transitions 
the CO emission reaches a minimum at the center. This is 
also seen in the higher-resolution array maps published by 
Casoli & Gerin (1993) and Garci'a-Burillo et al. (2003). 

In Fig.m the J=3-2 ^^CO distribution appears more 
compact than that in the J=2-l transition. This is partly 
an effect caused by the significantly higher resolution (14" 
vs 21") of the former. However, closer analysis shows that 
the higher-excited CO gas is indeed more concentrated than 
the CO in the ground states. In all observed transitions, CO 
emission extends out to i? = 30" (see e.g. Fig.[3]-top), but 
the half-width (FWHM) of the CO distribution, corrected 
for finite beam width, is 43" in the J=2-l transition, and 
only 20" (i? = 200 pc) in both the J=3-2 and J=4-3 tran- 
sitions (c/Fig.[9]-bottom), as well as in [CI] (not shown). It 
thus appears that the hotter molecular gas is preferentially 
found closer to the starburst center. 

The i^CO and ^^CO line ratios observed in NGC 4826 
are not very well fitted by the model parameters given in 
Table[B] In particular the J=4-3 ^^CO line is too strong, a 
situation we also found in the center of NGC 1068. The high 
temperature of this relatively dense (nH2 ~ 3000 cm~^ 
gas close to the nucleus opens up the possibility that X-ray 
emission from the latter plays a role in the excitation. A mi- 
nor radio source (Hummel et al. 1987; Turner & Ho 1994) 
in the center suggests the presence of a weak AGN (see also 
Garcfa-Burillo et al. 2003; Haan et al. 2008). It is notewor- 
thy that the molecular gas constitution changes away from 
the central peak. For instance, at a position 18" (355 pc) 
ESE from the nucleus we find line ratios that suggest ei- 
ther a kinetic temperature of about 30 K and a density 
of about 3000 cm""^, or comparable amounts of cold, very 
dense {Tkin = 10 K, = 10^ cm~^) and warm, very 

tenuous {Tkin = 60 K, riHa = 10^ cm~^) molecular gas. 
The lack of more ^'^CO measurements makes it impossi- 
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Fig. 9. Top: major-axis-velocity diagram of the J=3-2 
i^CO emission from NGC 4826. Contour step is 40 mK. 
Bottom: J=4-3 ^^CO map of NGC 4826 with contour step 
of 10 Kkms-i. 



ble to choose between these alternatives, but especially the 
first one is compatible with molecular gas heated by UV 
radiation from luminous stars. 



5.5. NGC 7469 

As expected, our maps (Fig. [51) do not show much detail 
in this distant galaxy. Nevertheless, the J=3-2 ^^CO maps 
show emission from the center superposed on a more ex- 
tended disk. Unfortunately, our J=4-3 ^^CO map is in- 
complete in the direction of the J=3-2 extension. Our an- 
gular resolution is insufficient to resolve the rotation curve 
of NGC 7469, but once again it can be clearly seen that the 
nucleus itself is at a relative CO minimum, as confirmed by 
the array maps obtained by Davies et al. (2004). We may 
thus expect the molecular gas parameters to be dominated 
by the extended starburst region rather than by the imme- 
diate surroundings of the AGN. The observed fine ratios are 
well-fitted by the model parameters listed in Tabled] These 
imply that a small amount (typically a few times 10^ Mq) of 
the molecular gas is modestly warm {Tkin = 30-60 K) and 
rather dense (nHa = 10^ cm~'^). It is tempting to asso- 
ciate this relatively small amount with the AGN-irradiated 
gas also seen in near-infrared H2 emission (Genzel et al. 
1995; Davies et al. 2004). The bulk of the molecular gas 
has mean temperatures {Tkin = 1-30 K and mean densi- 
ties (nH2 = 10^ — 10'^ cm"-^) appropriate to extended star 
formation in the disk and spiral arms. 
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5.6. Column densities and masses 

Characteristically, we find H2 column densities averaged 
over a 21" beam of roughly 1-3 x 10^^ cm^^, correspond- 
ing to total hydrogen columns of 2-6 x 10^^ cm~^. These 
total columns also include the ambient contribution by 
neutral atomic hydrogen which we have so far neglected. 
HI absorption column densities actually measured towards 
these galaxies often exceed the H2 column densities listed 
in TableH Towards NGC 1068 HI column densities N{HI) 
= (1-2) xlO^i cm"2 (Brinks et al. 1997) are seen, but to- 
wards NGC 2146 peaks are up to six times higher than our 
average H2 value, and reach as high as N{HI) — (5-12) 
xlO^i cm-2 (Tarchi et al. 2002) in a 2.1" x 1.5" beam. 
Towards the nucleus of the edge-on galaxy NGC 3079, 
Pedlar et al. (1996) even find N{HI) = 27 x lO^^ cm^^ 
in a 4" x 3" beam. Compared to this, the mostly face-on 
galaxy NGC 7469 has rather moderate neutral atomic hy- 
drogen column densities N{HI) « 3 x 10^^ cm~^), seen 
directly in absorption against the nucleus (Beswick et al. 
2002) and also implied by X-ray absorption (Barr, 1986). 
Especially in edge-on galaxies, not all HI seen in absorption 
necessarily resides in the central region we have sampled. 
Moreover, in all these cases the effective absorbing area 
is much smaller (by factors of 30 to 150) than the large 
21" beam over which we have averaged our H2 in Table[71 
Our models imply CO, hence H2, surface filling factors of 
5-8 for NGC 2146, NGC 4826, and NGC 7469, of 15 for 
NGC 1068, and 33 for NGC 3079. This suggests that the 
HI contribution to the total hydrogen column densities im- 
plied by the models is very small for NGC 1068 (on the 
order of a few per cent), and somewhat more (10-30 per 
cent) for NGC 2146, NGC 3079 and NGC 7469. In view of 
the uncertainty of a factor of two inherent in the model val- 
ues, the error introduced by neglecting the HI contribution 
is negligible. 

Our H2 masses are much lower than those calcu- 
lated previously by authors who used a 'standard' CO- 
to-H2 conversion factor X. For instance, in the center of 
NGC 2146 we find a molecular hydrogen mass ~ 
2 X 10^ M0, a factor of 23 lower than the mass listed 
by Greve et al. (2006). Even taking into account that 
these authors used the rather high conversion {X = 3 x 
10^° cm~^/( K kms~^)), there is still an order of magnitude 
difference with our result. In the galaxy centers discussed 
here, we find X = 0.1-0.2 x 10^° cm-V(Kkms"i) which 
is a factor of 10-20 lower than the 'standard' Milky Way 
X factors (see also Strong et al. 2004), but quite similar 
to those found in other galaxy centers (see Israel & Baas 
2003 and references therein). Clearly, X factors derived for 
molecular clouds in the Solar Neighborhood do not apply 
to molecular gas in the very different environment posed 
by galaxy centers. We also note that the present results are 
consistent with the low values derived independently for 
the same galaxy centers by other authors from e.g. analysis 
of infrared observations (c/Genzel et al. 1995; Braine et al. 
1997; Papadopoulos & Allen, 2000; Usero et al. 2004; Davies 
et al. 2004). Using the 'standard' value of X, authors in the 
past have concluded to very high gas mass fractions (up to 
half of the dynamical mass) in the centers of galaxies. With 
the revised low values, there is no longer a need to assume 
such extreme conditions. In fact, in our galaxies we find 
much more modest mass ratios Mn^/Mdyn ~ 0.5 — 1.0%. 



5. 7. Temperature and excitation 

Although a classical UV-illuminated photon-dominated re- 
gion (PDR) associated with a star-forming region may 
reach quite high kinetic temperatures, the amount of gas 
at such elevated temperatures is rather small and limited 
to the immediate vicinity of the stellar heating source. It is 
very difficult, even in intense star-forming regions, to reach 
globally averaged temperatures in excess of 20-30 K. An 
example will illustrate this. If the heating is provided by 
luminous (L = 0.5 - 2.6 x 10*^ Lq. Vacca et al. 1996) 05 
stars, and perfect temperature equilibrium is assumed, the 
separation between star and heated gas should be about 
1.5-2.5 pc at 20 K, and only 0.6-1.0 pc at 30 K. 

When we look at TablelHlwe see that only the molecu- 
lar line emission from the clear Seyfert galaxies NGC 1068 
and NGC 7469 may arise from molecular gas dominated by 
a starburst. In the previous section, we have seen that in 
either galaxy our resolution was insufficient to separate the 
emission from the circumnuclear molecular gas from that of 
the much more extended molecular gas associated with the 
more starburst. We may therefore expect that observations 
similar to those presented in this paper, but conducted at 
substantially higher resolutions, might indeed require a dif- 
ferent excitation for the gas close to the nucleus. As also 
noted above, there are already indications that this so in 
NGC 1068. In this respect, millimeter-array observations of 
this galaxy in high-J transitions of both ^^CO and ^'^CO 
(SMA, ALMA) would be quite interesting. 

Much if not all of the molecular gas in NGC 2146, which 
appears to be a starburst- dominated galaxy, is in fact far 
too hot {Tkin = 125-150 K) to be excited by UV photons 
from luminous stars. The situation is more ambiguous in 
the merger galaxy NGC 4826, where two thirds of the mass 
is probably heated in PDR fashion, but one third of the gas 
is again too hot with Tkin ~ 60 K. In NGC 3079, 20 - 60% 
of the gas may be heated by the starburst that is present at 
a few hundred parsec from the nucleus, but the remaining 
relatively diffuse gas is far too hot with a temperature as 
high as Tkin ~ 150 K. It is noteworthy that the two galax- 
ies with the hottest gas (NGC 2146 and NGC 3079) are 
both associated with a wind-driven X-ray outflow (Armus 
et al. 1995; Delia Ceca et al. 1999; Pietsch et al. 1998). 
One might speculate that the hot molecular gas compo- 
nent is efficiently heated by the X-rays from these flows 
(see Meijerink et al. 2006, 2007) rather than by the UV 
photons from the starburst. 

6. Conclusions 

1. We have observed the centers of five galaxies with sig- 
nificant central activity in the first four transitions of 
^^CO and the first three transitions of ^'^CO, and three 
of them also in the lower transition of [CI] . 

2. All galaxies were mapped over at least 1' x 1' in the 
J=2-l and J=3-2 and over smaller areas in the J=4-3 
transitions of ^^CO. 

3. The line ratios, reduced to a standard 21" beam vary 
significantly in the sample galaxies. In general, the J—A- 
3 line is still relatively strong, indicating elevated tem- 
peratures for the molecular gas. 

4. Optical depths in ^^CO appear not to be very high, 
with ^^CO/^'^CO ratios generally in the range of 11- 
16. Exceptions are NGC 4826 and NGC 7469. The first. 
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a merger galaxy, has isotopic ratios in the range 6-8 sug- 
gesting somewhat higher optical depths reminiscent of 
extended molecular clouds in the Solar Neighborhood. 
On the other hand, the distant Scyfcrt galaxy NGC 7469 
has isotopic ratios in the range of 16-20, suggesting that 
much of the gas in this galaxy has a low optical depth 
characteristic of relatively tenuous gas. 

5. The molecular gas is present in the form of clearly 
identifiable central condensations of outer radius R = 
0.3 — 1.5 kpc. In all galaxies, the molecular gas distri- 
bution exhibits a local minimum at the nucleus itself. 

6. In all galaxy central regions, we find moderate H2 
masses ofl-SxlO^Mo, except in the nearby NGC 4826 
where the mass is no more than 3 x IO'^'Mq. These 
masses are typically about one per cent of the dynamical 
mass in the same region. The CO-to- H2 conversion fac- 
tors X implied by these masses and the observed ,7=1-0 
"'^^CO line intensities is typically an order of magnitude 
lower than the 'standard' Solar Neighborhood X-factor. 

7. Intriguingly, the excitation of molecular gas in the 
Seyfert galaxies NGC 1068 and NGC 7469 appears to be 
dominated by UV-photons from the extended circum- 
nuclear starburst. In contrast, the molecular gas in the 
starburst galaxy NGC 2146, and part of the molecular 
gas in NGC 3079 (starburst cum AGN) are probably ex- 
cited by Xrays from their massive wind-driven outflows. 
The situation for NGC 4826 is less clear. 

Acknowledgements. I am much indebted to the late Pred Baas, who 
died suddenly and far too early in 2001, for carefully planning and 
performing many of the observations described in this paper. I also 
thank the various JCMT observers who made additional observations 
in queue mode. Finally, it is always a pleasure to thank the personnel 
of both the JCMT and the IRAM 30m telescope for their able and 
generous support. 

References 

Armus L., Shupe D.L., Matthews K., et al. 1994, in: 'Infrared 
Astronomy with Arrays', ed. I. McLean, Kluwer (Dordrecht), p. 
501 

Armus L., Heckman T.M., Weaver K.A., & Lehnert M.D., 1995, ApJ 
445, 666 

Bayet E., Gerin M., Phillips T.G., & Contursi A., 2004, A&A 427, 45 
Beswick R. J., Pedlar A., & McDonald A.R., 2002 MNRAS 335, 1091 
Blustin A.J., Kriss G.A., Holezer T., et al. 2007, A&A 466, 107 
Braun R., Walterbos R.A.M., Kennicutt R.C., & Tacconi L.J., 1994, 

ApJ 420, 558 

Brainc J., Guclin M., Dumkc M., ct al. 1997, A&A 326, 963 
Braun R., Walterbos R.A.M., & Kennicutt R.C., 1992, Nature 360, 
442 

Braun R., Oostcrloo T.A., Morganti R., Klein U., & Beck R., 2007, 
A&A 461, 455 

Blietz M., Cameron S., Drapatz S., et al. 1994 ApJ 421, 92 

Brinks E., Skillman E. D., Terlevich R. J., & Terlevich E., 1997 

Ap&SS248, 23 
Casoli P., & Gerin, 1993, A&A 279, L41 

Cecil C, Blad-Hawthorn J., Veilleux S.,& Pilippcnko A.V., 2001 ApJ 
555, 338 

Condon J. J., Condon M.A., Gislcr G., & PuschcU J. J., 1982, ApJ 252, 
102 

Condon J. J., Hclou G., Sanders D.B., & Soifcr B.T., 1990 ApJS 73, 
359 

Condon J. J., Hclou G., Sanders D.B., & Soifcr B.T., 1996 ApJS 103, 
81 

Crawford M.K., Gcnzcl R., Towncs C.H. & Watson D.M., 1985 A&A 
291, 755 

Davics R.I., Sugai H., & Ward M.J., 1998 MNRAS 295, 43 
Davics R.I., Tacconi L.J., & Genzel R., 2004 ApJ 602, 148 
de Bruyn A.G., 1977, A&A 58, 221 

Dumke M., Nieten Ch., Thuma G., Wielebinski R. & Walsh W., 2001 
A&A 373, 853 



Delia Ccca R., Griffiths R.E., Heckman T.M., Lehnert M.D., Weaver 

K.A., 1999, ApJ 514, 772 
Fabbiano G., Kim D.-W., & Trinchicri G., 1992 ApJS 80, 531 
Gallimorc J. P., Baiim S.A., O'Dca CP., Brinks E., & Pedlar A., 1994, 

ApJL 422, L13 

Garcfa-Burillo S., Combes P., Hunt L.K., ct al., 2003 A&A 407, 485 
Garnctt D.R., Shields G.A., Pcimbcrt M., et al. 1999 ApJ 513, 168 
Genzel R., Weitzcl L., Tacconi-Garman E., ct al. 1995, ApJ 444, 129 
Gerin M., & Phillips T.G.., 2000, ApJ 537, 644 

Grcvc A., Nciningcr N., Sievers A., & Tarchi A., 2006, A&A 459, 441 
Haan S., Schinncrcr E., Mundell C.G., Garcia-Burillo S., & Combes 

P., 2008 AJ 135, 232 
Hawardcn T.G., Israel F.P., Geballe T.R., & Wade R., 1995, MNRAS 

376, 1197 

Hclfcr T.T., & Blitz L., 1995 ApJ 450, 90 

Hclfcr T.T., Thornlcy M.D., Regan M.W., ct al. 2003, ApJS 145, 259 
Hcnkcl C, Maucrsbcrgcr R., Wiklind T., ct al. 1993, A&A 268, L17 
Henkel C, Wliiteoak J.B., & Mauersberger R., 1994 A&A 284, 17 
Henkel C, Chin Y.-N, Mauersberger R. & Whiteoak J.B., 1998 A&A 

329, 443 

Hummel E., van der Hulst J.M., Keel W.C., & Kennicutt R.C., 1987 

A&AS 70, 517 

Hogcrhcijdc M.R., & van der Tak F.P.S., 2000 A&A 362, 697 
ImpcUizzeri C.M.V., Hcnkcl C., Roy, A.L., & Mcnten K.M., 2008, 
A&A 484, L43 

Inui T., Matsumoto H., Tsuru, T.G., Koyama K., Matsushita S., & 

Peck A.B., 2005 PASJ 57, 135 
Irwin J. A., & Saikia D.J., 2003, MNRAS 344, 977 
Israel P.P., White G.J., & Baas P., 1995, A&A 302, 343 
Israel P.P., van der Wcrf P.P., Hawardcn T.G., & Aspin C, 1998 A&A 

336, 433 

Israel P.P., & Baas P., 1999, A&A 351, 10 
Israel P.P., & Baas P., 2001, A&A 371, 433 (Paper I) 
Israel P.P., & Baas P., 2002, A&A 383, 82 
Israel P.P., & Baas P., 2003, A&A 404, 495 (Paper II) 
Israel P.P., Tilanus R.P.J. , & Baas P., 2006, A&A 445, 907 (Paper 
III) 

Irwin J. A., & Scaquist E. R., 1991 ApJ 371, 111 

Irwin J. A., & Sofuc Y., 19921 ApJL 396, L75 

Jackson J.M., & Ho P.T.P., 1988, ApJL 324, L5 

Jackson J.M., Paglione T.A.D., Ishizuki S., & Nquyen-Q-Rieu, 1993, 

ApJL 418, L13 
Janscn D.J., 1995, Ph.D. thesis. University of Leiden (NL) 
Janscn D.J., van Dishoeck E.F. & Black J.H., 1994, A&A , 282, 605 
Kaneko N., Morita K., Pukui Y., and 5 other authors, 1989, ApJ 337, 

691 

Kaneko N., Morita K., Pukui Y., and 4 other authors, 1992, PASJ 44, 
341 

Koda J., Sofuc Y., Kolmo K., ct al., 2002, ApJ 573, 105 
Kronberg P.P., & Bicrmann, P., 1981, ApJ 243, 89 
Le Ploc'h E., Mirabel L.P., Laurent O., et al. 2001 A&A 367, 487 
Lepp S., & Dalgarno A., 1996 A&A , 306, L21 

Liscnfcld U., Alexander P., Pooley6 G.G., & Wilding, T. 1996 MNRAS 

281, 301 

Lumsdcn S.L., Moore, T.J.T., Smith C, ct al. 1999 MNRAS 303, 209 
Maucrsbcrgcr R. & Hcnkcl C., 1993 Rev. Mod. Astron. 6, 69 
Meijcrink R., Spaans M., & Israel P.P., 2006 ApJL 650, L103 
Mcijcrink R., Spaans M., & Israel P.P., 2007 A&A 461, 793 
Mcixncr M., Puchalsky R., Blitz L., Wright M., & Heckman T., 1990 
ApJ 354, 158 

Moustakas J., & Kennicutt R.C., 2006 ApJ 651, 155 
Papadopoulos P.P., & Scaquist E.R., 1999, ApJ 516, 114 
Paglione T.A.D., Wall W. P., Young J.S. ct al. 2001, ApJ 135, 183 
Pedlar A., Mundell C.G., Gallimorc J. P., Baum S.A., & O'Dea CP., 

1996, VA, 40, 91 
Pictsch W., Trinchicri C, & Vogler A., 1998, A&A 340, 351 
Pilyugin L.S., Vi'lchez J.M., & Contini T., 2003 A&A 425, 849 
Planesas P., Gomez-Gonzalez J., & Martin-Pintado J., 1989, A&A 

216, 1 

Planesas P., Scoville N.Z., & Myers S.T., 1991, ApJ 369, 364 

Rix H.-W., Kennicutt R.C, Braun R., & Walterbos R.A.M., 1995, 

ApJ 438, 155 
Rosolowsky E., & Blitz L., 2005, ApJ 623, 826 
Rotaciuc v., Krabbc A., Cameron M., et al. 1991 ApJL 370, L23 
Rubin V.C, 1994, AJ 107, 173 

Sakamoto K., Okumura S.K., Ishizuki S. & Scoville N.Z., 1999 ApJS 
124, 403 



F.P. Israel: CI and CO in five bright spirals 



15 



Sandagc A. & Bcdkc J., 1994, The Carnegie Atlas of Galaxies. Volume 

I, Carnegie Institution of Washington 
Schinnerer E., Eckart A., Tacconi L.J., Genzel R., & Downes D., 2000, 

ApJ 533, 850 

Scoville N.Z., Young J. S. & Lucy L.B., 1983, ApJ 270, 443 
Sil'chenko O.K., 1996, AL 22, 108 
Sofue Y., & Irwin J.A., 1992, PASJ 44, 353 
Sofue Y., Koda J., Kohno K., et al. 2001, ApJL 547, L115 
Sternberg A., Genzel R., & Tacconi L., 1994, ApJL 436, L131 
Stacey G.J., Geis N., Genzel R., et al. 1985 A&A 373, 423 
Strong A.W., Moskalenko LV., Reimer O., Digel S., & Diehl R., 2004 
A&A 422, 47 

Tacconi L.J., Genzel R., Blietz M., Cameron M., Harris A.L, & 

Madden S., 1994, ApJL 426, L77 
Tacconi L.J., BUctz M., Cameron M., et al., 1996, VA 40, 23 
Tarchi A., Neininger N., Grcvc A., et al. 2000, A&A 358, 95 
Tarchi A., Henkcl C, Peck A.B., & Mentcn K.M., 2002, A&A 389, 39 
Tarchi A., Greve A., Peck A.B., et al. 2004, MNRAS 351, 339 
Tinncy C.G., Scoville N.Z., & Soifcr B.T., 1990, ApJ 362, 473 
Trotter A.S., Grecnhill L.J., Moran J.M., c3t al. 1998, ApJ 495, 740 
Turner J.L. & Ho P.T.P., 1994 ApJ 421, 122 

Usero A., Garcfa-Burillo, Fuentc A., Marti'n-Pintado J., & Rodrii'guez- 

Fernaandez N.J., 2004, A&A 419, 897 
Vacca W.D., Garrnany CD., & ShuU J.M., 1996, ApJ 460, 914 
van Dishocck E.F. & Black J.H., 1988, ApJ 334, 771 
Vcilleux S., Cecil G., Bland-Hawthorn J., et al. 1994, ApJ 433, 48 
Vila-Costas M.B. & Edmunds M.G. 1992 MNRAS 259, 121 
Walterbos R.A.M., Braun R.,& Kennicutt R.C., 1994, AJ 107, 184 
Wilson A.S., & Ulvestad J.S., 1987 ApJ 319, 105 
Wilson A.S., Baldwin J. A., Sun S.-D., & Wright A.E., 1986, ApJ 310, 

121 

Wilson A.S., Helfcr T.T., Haniff C.A., & Ward M.J., 1991, ApJ 381, 
79 

Wynn-WiUiams C.G., BeckUn E.E., & Scoville N.Z., 1985, ApJ 297, 
607 

Xie S., Young J.S., & Schlocrb F.P., 1994, ApJ 421, 434 
Young J.S., Claussen M.J., & Scoville N., 1988a ApJ 324,115 
Young J.S., Claussen M.J., Kleimnaim S.G., Rubin V.C., & Scoville 

N.Z., 1988b ApJL 331, L81 
Young J.S., Wilson A.S., & ShopbcU P.L., 2001 ApJ 556, 56 
Young J.S., Xie S., Tacconi L., et al. 1995 ApJS 98, 219 
Zaritsky D., Kennicutt R.C. & Huchra J. P., 1994, ApJ 420, 87 
Zhao J.-H., Anantharamaiah K.R., Goss W.M., & Viallefond F., 1996, 

ApJ 472, 54 



